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Abstract. We study the coarsening dynamics of the three-dimensional random field 
Ising model using Monte Carlo numerical simulations. We test the dynamic scaling and 
super-scaling properties of global and local two-time observables. We treat in parallel 
the three-dimensional Edward- Anderson spin-glass and we recall results on Lennard- 
Jones mixtures and colloidal suspensions to highlight the common and different out of 
equilibrium properties of these glassy systems. 

1. Introduction 

The physics of domain growth is well understood [U Just after the initial thermal 
quench into the ordered phase, the spins in a ferromagnetic system tend to order 
and form domains of the equilibrium states. In clean systems the ordering dynamics 
is governed by the symmetry and conservation properties of the order parameter. 
When impurities are present the dynamics are naturally slowed down by domain-wall 
pinning [3j. The dynamic scaling hypothesis states that the time-dependence in any 
macroscopic observable enters only through a growing length scale, R(t), either the 
instantaneous averaged or typical domain radius. However, a complete description of the 
phenomenon is lacking. In the pure cases the scaling functions are not known analytically 
and no fully satisfactory approximation scheme to estimate them is known [lj. In 
presence of disorder the limitations are more severe in the sense that the growth laws 
are derived by assuming that the relaxation is driven by activation over free-energy 
barriers and the properties of the latter are estimated with energy balancing arguments 
applied to single interfaces that are hard to put to the test. Even in the relatively simple 
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random bond Ising model (RBIM) the time dependence of the growth law remains a 
subject of controversy [I]. 

Quenched randomness may be weak or strong in the sense that the first type does 
not change the nature of the low-temperature phase, as in the random bond or random 
field Ising (RFIM) models, or it can change it as in spin-glasses. Fisher and Huse 
conjectured that in the first class of systems, once the scaling hypothesis is used to 
describe the long times dynamics, so that times and lengths are measured in units of 
R(t), no out of equilibrium observable depends on the quenched randomness [5] and 
their scaling functions are thus identical to the ones of the pure limit. This is the so- 
called 'super-universality' hypothesis in coarsening phenomena. Tests of this hypothesis 
as applied to the equal-times two-point function of the 3d RFIM and the 2d RBIM 
appeared in [6] and [Tj, respectively, and the distribution of domain areas in the 2d 
RBIM in |g]. 

The dynamics of generic glassy systems is less well understood but presents some 
similar aspects to those mentioned above. The droplet model of finite-dimensional spin- 
glasses is based on the assumption that in the low-temperature phase these systems 
also undergo domain growth of two competing equilibrium states [5] . In the mean-field 
limit spin-glasses have, though, a very different kind of dynamics [9j [10] that cannot 
be associated to a simple growth of two types of domains. Numerical studies of the 3d 
Edwards- Anderson (EA) model [HJ [121 [321 El EE] have not been conclusive in deciding 
for one or the other type of evolution and, in a sense, show aspects of both. A one- 
time dependent 'coherence '-length, R(t), has been extracted from the distance and time 
dependence of the equal-time overlap between two replicas evolving independently with 
the same quenched disordered interactions [321 HH [15]. A power-law R(t) ~ t 1 l z ^ r ) 
with the dynamic exponent z(t) = z(T c )T c /T fits the available data for the 3d EA 
and z(T c ) = 6.86(16) with Gaussian [15] and z(T c ) = 6.54(20) with bimodal [121 02] 
couplings. Still, it was claimed in [15] that the overlap decays to zero as a power law at 
long distances and long times such that r/R(t) is fixed, implying that there are more 
than two types of growing domains in the low temperature phase. 

A two-time dependent length, t w ), can be extracted from the analysis of the 
spatial decay of the correlation between two spins in the same system at distance r 
and different times t and t w after preparation [16] . The latter method is somehow 
more powerful than the former one in the sense that it can be easily applied to glassy 
problems without quenched disorder. If there is only one characteristic length-scale in 
the dynamics R(t) should be recovered as a limit of £,(t,t w ) but this fact has not been 
demonstrated. 

The mechanism leading to the slow relaxation of structural glasses is also not 
understood. Still, molecular dynamic studies of Lennard- Jones mixtures [TF] and the 
analysis of confocal microscopy data in colloidal suspensions [18] show that two-time 
observables have similar time dependence as in the 3d EA model. Two-time correlations 
scale using ratios of one-time growing functions that, however, cannot be associated to 
a domain radius yet. A two-time correlation length £ with characteristics similar to the 
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one in the 3d EA can also be denned and measured. 

The understanding of dynamic fluctuations in out of equilibrium relaxing systems 
appears as a clear challenge [I9j. In systems with quenched randomness different sample 
regions feel a different environment and one expects to see their effect manifest in 
different ways working at fixed randomness. The effect of quenched randomness is 
at the root of Griffiths singularities in the statics of disordered systems, for instance. In 
structural or polymer glasses there are no quenched interactions instead, but still one 
expects to see important fluctuations in their dynamic behaviour both in metastable 
equilibrium and in the glassy low temperature regime. The question of whether the 
fluctuations in generic glassy systems resemble those in coarsening systems has only 
been studied in a few solvable cases such as the model of ferromagnetic coarsening in 
the large N limit [20] and the Ising chain [21] . 

In this paper we study ferromagnetic ordering in the 3d RFIM following a quench 
from infinite temperature and we compare it to the dynamics of the 3d EA spin- 
glass and particle glassy systems. Our aim is to signal which aspects of their out 
of equilibrium evolution differ and which are similar by focusing on freely relaxing 
observables - no external perturbation is applied to measure linear responses. We 
test the scaling and super-universality hypothesis in the RFIM and we explicitly show 
that the latter does not apply to the EA model. We analyse the spatio-temporal 
fluctuations in the coarsening problem and we compare them to the ones found in 
spin-glasses [T^l [16] . the O(N) ferromagnetic coarsening in the large N limit [20], and 
other glassy systems [T8(, [22l [23] . 

The organisation of the paper is the following. In Sect. [2] we define the models and 
we describe the numerical procedure. Section [3] is devoted to the study of the growing 
length scale, R, the scaling and super-universality hypothesis, and the two-time growing 
length, £. In Sect. 0]we focus on the local fluctuations of two time observables. We study 
two-time coarse-grained correlations and we analyse their statistical properties as time 
evolves. Finally, in Sect. [5] we present our conclusions. 

2. The models 

Two varieties of quenched disorder are encountered in spin models: randomness in the 
strength of an externally applied magnetic field and randomness in the strength of the 
bonds. The RFIM and the EA spin-glass are two archetypal examples of these. In this 
Section we present their definitions and we recall some of their main properties. 

2.1. The Random Field Ising Model 

The 3d Random Field Ising model (RFIM) is defined by the Hamiltonian [21] 




(1) 
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The first term encodes short range ferromagnetic (J > 0) interactions between nearest 
neighbour Ising spins, = ±1, placed on the nodes of a cubic lattice with linear size L. 
Hi represents a local random magnetic field on site %. We adopt a bimodal distribution 
for these independent identically distributed random variables [Hi = ±H with equal 
probability). H quantifies the strength of the quenched disorder. Hereafter we set J = 1 
and we use units in which k B = 1. 

The RFIM is relevant to a large class of materials due to the presence of defects 
that cause random fields. Dilute anisotropic antiferromagnets in a uniform field are the 
most studied systems expected to be described by the RFIM. Several review articles 
describe its static and dynamic behaviour [3] and the experimental measurements in 
random field samples have been summarized in [25]. Dipolar glasses also show aspects 
of random field systems [26] . 

In the case H = 0, the RFIM reduces to the well known clean Ising model with a 
phase transition from a paramagnetic to a ferromagnetic state occurring at T c ~ 4.515. 
It is well established that in d = 3 (not in d — 2) there is a phase separating line on the 
(T, H) plane joining (T c , H = 0) and (T = 0, H c ). At T = and small magnetic field, 
it has been rigorously proven that the state is ferromagnetic [271 -28J . The nature of the 
transition close to zero temperature has been the subject of some debate. Claims of it 
being first order [29] have now been falsified and a second order phase transition has 
been proven [301 El] • It was also argued but not established that there might be a spin- 
glass phase close to (T = 0, H c ) [32]. Recent numerical simulations yield H c ~ 2.215(35) 
at T = [331 [31]. 

2.2. The Edwards- Anderson spin-glass 

The 3d Edwards-Anderson (EA) spin glass is defined by 



The interaction strengths act on nearest neighbours on a cubic three-dimensional 
lattice and are independent identically distributed random variables. We adopt a 
bimodal distribution, J^- = ±1 with equal probability. This model undergoes a static 
phase transition from a paramagnetic to a spin-glass phase at T g ~ 1.14(1) [35J. 
The nature of the low temperature static phase is not clear yet and, as for the out 
of equilibrium relaxation, two pictures developed around a situation with only two 
equilibrium states as proposed in the droplet model and a much more complicated 
vision emerging from the solution of the Sherrington-Kirkpatrick model, its mean-field 
extension [36] . 

2.3. Numerical methods 




(2) 



We focus here on the out of equilibrium relaxation in the ordered phase. We simulate 
the dynamics following an instantaneous quench from infinite temperature at the initial 
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time, t = 0, by choosing a random initial condition: Sj(t = 0) = ±1 with probability one 
half. The order parameter is not conserved during the evolution. We use the continuous 
time Monte Carlo (MC) procedure [371 1381 139]. This algorithm, which is nothing else 
than a re-organisation of the standard Monte Carlo rule, is rejection free. This makes it 
spectacularly faster than standard MC which would have a rejection rate close to 1 in 
the ferromagnetic phase of the RFIM. Times are expressed in usual Monte Carlo steps 
(MCs): 1 MCs corresponds to N = L d spin updates with the standard MC algorithm. 
The way to translate from the continuous time MC to standard MC units, in which we 
present our results, is explained in [37J [381 EH] ■ 

We study the relaxation dynamics with non-conserved order parameter in the 
(d = 3) ferromagnetic phase of the RFIM at relatively low temperature and small 
applied field. Interesting times are not too short - to avoid a short transient regime - 
and not too long - to avoid reaching equilibration (in ferromagnetic coarsening a non- 
zero magnetization density indicates that the coarsening regime is finished and other 
more refined methods are used in the spin-glass case [10]). We delay equilibration by 
taking large systems since the equilibration time rapidly grows with the size of the lattice. 
A reasonable numerical time-window is [10 3 , 10 7 ] MCs. We show results obtained using 
lattices with L = 250 (N = 1.5 x 10 7 spins) in the RFIM and L = 100 (N = 10 6 spins) 
in the spin-glass. We checked that finite size effects are not important in any of these 
cases for averaged quantities. 

3. The typical growing length 

In this Section we study the typical growing length (a geometric object) in the RFIM 
and the EA model. We establish scaling and super-universality relations for three types 
of correlations functions (statistical objects). Two of them involve either two space 
points and one time, or one space point and two times, and are the usual observables 
studied in coarsening phenomena. The third one is commonly used in the study of glassy 
systems where two-point correlations are not sufficient to characterize the dynamics of 
the systems [HI [HU [ITJ [18] and allows for the definition of a two-time dependent length 
that we can compare to the one obtained in the 3d EA model and glassy particle systems. 

3.1. The RFIM 

During the ferromagnetic coarsening regime, there are as many positive as negative spins 
in such a way that the magnetization density stays zero in the thermodynamic limit and 
weakly fluctuates around zero for finite size systems. Everywhere in the sample, there 
is a local competition between growing domains. Eventually, after an equilibration time 
r eq (that diverges with the system size), one of the two phases conquers the whole system 
scale. 

In the coarsening regime (times shorter than r eq ) dynamic scaling [1] applies and the 
growth of order is characterized by a typical domain radius, R(t;T,H), that increases 
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in time and depends on the control parameters, T and H, and the dimension of space, 
d [H]. While in the absence of impurities it is clearly established that, for non-conserved 
order parameter dynamics, the domain length R grows as R ~ t}l 2 independently of 
d [1] with a prefactor that monotonically decreases upon increasing temperature |42j, 
the functional form of R is less clear in random cases. Scaling arguments based on the 
energetics of single interfaces j3l H3j HH US] predict a crossover from the pure case result 
at short time-scales when it is easy to inflate, to a logarithmic growth, 

R(t;H,T) = ^\n(t/r(T,H)) . (3) 

The fact that the prefactor grows with T (as opposed to what happens for pure curvature 
driven dynamics [12]) is due to the activated character of the dynamics. Several 
proposals for the characteristic time r exist: r ~ (T/H 2 ) 2 [4~3l 12] and r ~ Toe A<yT ^ H2 
with A(T) a weakly temperature dependent function [Bj. To ease the notation in what 
follows we do not write explicitly the T and H dependence of R. 

From the point of view of the renormalization group, all points within the 
ferromagnetic region of the (T, H) phase diagram flow to the stable, zero-temperature, 
zero-disorder sink. Hence, randomness and temperature should be irrelevant in 
equilibrium at T < T c . The super-universality hypothesis states that for non equilibrium 
ordering dynamics, once lengths are scaled with the typical length R, quenched random 
fields are irrelevant and all scaling functions are the ones of the pure 3d Ising system at 
T = with non-conserved order parameter. 

3.1.1. The equal-time spatial correlation. A careful analysis of the field and time 
dependence of the growing length scale together with tests of the scaling hypothesis 
applied to the equal-time correlation 

C 2 (r;t) = ( Sl (t) Sj (t))^_ rll=r , (4) 

where the average runs over all spins in the sample, appeared in [6JHH]. In the coarsening 
regime, at distances a <^ r <^ L with a the lattice spacing and r/R(t) finite, C 2 (r; t) is 
expected to depend on r and time t only through the ratio r/R, 

C 2 (r;t)^m 2 eq f 2 (r/R(t)), (5) 

with m eq the equilibrium magnetization density (that decreases with increasing T and/or 
H), limjc-^o — 1 an d lim^oo f 2 (x) = 0. Since the spatial decay is approximately 
exponential, C 2 (r; t) oc e~ r ^ R ^ for not too long r, we use this functional form to extract 
R from the data fit at each set of parameters (T,H,t). Figured] (a) shows that the 
growing length R has two regimes: shortly after the quench R grows as t 1 ^ 2 like in the 
pure case and it later crosses over to a logarithmic growth. This is consistent with 
previous numerical studies in 2d (7J H7] and 3d systems (6J 08]. In Fig. [1] (b) we test the 
dependence on T and H by plotting ^R versus t/r for T = 1, 2 and H = 0.5, 1, 1.5. 
We found the best collapse using r ~ but the precision of our data is not high 
enough to distinguish between this and the rs proposed in [13] and [6j. Our numerical 
results tend to confirm the T/H 2 dependence of R even in the early stages of the growth. 
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Figure 1. (a) With line-points (red), the growing length R(t) at T — 1 and H = 1. 
The green curve is the power law \fi that describes well the data at short times, right 
after the temperature quench. The blue line is a logarithmic law apt to describe the 
behaviour at longer time-scales. In the inset: the same data in a log-log scale to 
highlight the quality of the \fi behaviour at short times, (b) Study of the dependence 
of R on the parameters T and H for two values of T and three random field strengths 
H given in the key. 




0.2 0.4 0.6 0.8 1 1.2 1.4 0.2 0.4 0.6 0.8 1 1.2 1.4 

r/R(t) r/R(t) 

Figure 2. (a) The scaling function f 2 (r/R) for T = 0.5, 1, 2 and H = 1, 1.5. (b) The 
same data in a linear- log scale showing that /2 is close to an exponential at short r/R. 



Since the work of [6], it is now clear that f% in Eq. (jHJ) is independent of H, and very 
similar to the one of the pure system. In Fig. [2] we also find that the scaling functions 
fi at different T fall on top of one another. Thus fi is independent of H and T. 

3.1.2. The two-time self-correlation. It is commonly defined as 

1 N 

C{t,t w ) = -J^M) 8 ^)) ^ (6) 

i=i 

and quantifies how two spin configurations of the same system, one taken at t w (waiting 
time) and the other one at t > t w , are close to each other. The angular brackets here 
indicate an average over different realizations of the thermal noise. In the large N limit, 
this quantity is self- averaging with respect to noise and disorder induced fluctuations. 
This two-time function has been used as a clock for the out of equilibrium dynamics of 
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glassy systems [HJ [TU] and we shall use this property again, in the study of the two-time 
growing length and fluctuations. 

The behaviour of C is well understood for coarsening systems. As long as the 
domain walls have not significantly moved between t w and t(> t w ) (that defines what 
we shall call later short time delay), the self-correlation is given by the fluctuations of 
spins that are in thermal equilibrium inside the domains. As any other equilibrium two- 
time function, the self-correlation depends then only on t — t w . Later, for longer time 
delays, the displacement of domain walls cannot be neglected any more and C looses its 
time-translational invariance. The self-correlation can be written as a sum of two terms 
representing the thermal and aging regimes: 

C{t, t w ) = C t h(t — t w ) + C ag (t, t w ) (7) 

with the limit conditions 

Cth(O) = 1 - <? EA , \im tw ^ t - C ag (t, t w ) = q EA , 

lim^^oo C th (t - t w ) = , lim t » iro C ag (t, t w ) = . 

Qea is a measure of the order parameter and in a ferromagnetic phase it simply equals 
m 2 eq , the magnetization squared. 

In Fig. [3] (a) we show the decay of the two-time correlation C as a function of the 
time delay t — t w for t w = 10 3 , 10 4 , 10 5 at T = 1 and H = 1. On each of these curves, 
one can distinguish the two dynamic regimes. The longer the waiting time the later the 
aging regime appears. In Fig. [3] (b) we show the decay of the two-time correlation as a 
function of time-delay for t w = 10 3 and five pairs of parameters (T, H) given in the key. 
It is clear that the full relaxation depends strongly on the external parameters: raising 
the temperature or reducing the random field strength speeds up the decay. For these 
values of T and H, q-^A does not change much but the decay in the aging regime does. 

Dynamic scaling implies that in the aging regime 

C ag (t,t w ) = q EA f ( -^ , (8) 

with R the typical length extracted from C2, /(l) = 1 and /(oo) = 0. For our choice of 
parameters (T, H), q E & is close to unity so we can easily compute / from the measured 
C by using / = C ag /qEA — C/<Zea- Super-universality states that / does not depend 
on T and H. In Fig. H] we show that both hypotheses apply to this quantity. In panel 
(a) we use a linear-linear scale while in panel (b) we present the same data in a double 
logarithmic scale. Although the scaling function / looks like a power law it is not. One 
expects that its tail [R(t) 3> R(t w )} becomes a power-law with an exponent A. The 
actual function / is not known. Most of the analytic efforts in domain growth studies 
are devoted to develop approximation schemes to derive /, / 2 and other scaling functions 
but none of them is fully successful pp. 



3.1.3. The four point- correlation function. In order to successfully identify a growing 
correlation length in glassy systems including the 3d EA spin-glass, one defines the 
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R(t)/R(U R(t)/R(t w ) 

Figure 4. Test of the scaling and super- universality hypothesis, (a) / = C ag /qEA vs. 
R(t)/R(t w ) at various pairs of (T,H) and t w given in the key. (b) The same data in 
log-log scale. 



two-time two-site correlation function [TU [16], [T71 HH 09] 

C 4 (r;t,t w ) = (si(t)s i (t w )s i (t)s i (t t0 ))| r - i _ r;j | =r . (9) 

We extract £ from its approximate spatial exponential decay: C^(r;t,t w ) — C 2 (t,t w ) oc 
e -r/i{t,t w ) relatively short r/^. (Other methods, like defining the connected four spin- 
correlation and extracting £ from its volume integral yield similar qualitative results 
though slightly different quantitatively.) Results of this analysis are shown in Fig. (a) 
where we plot £(t,t w ) as a function of t for different t w at T — 1 and H — 1. We 
identify a short t — t w regime that is independent of t w (thermal regime), whereas for 
long t — t w , time-translational invariance is broken (aging regime). In Fig. [5] (b) we plot 
€(t,t w ) versus 1 — C(t,t w ) for the three same values of t w , using t as a parameter. The 
dependence on 1 — C and t w is monotonic and very similar to the one obtained in the 
3d EA model [H] (see Fig. [7j). The thermal regime is almost invisible here since it is 
contained between C = 1 and C = gEA, with q^A — 1 for this set of parameters. We 
then propose 

£(t,t w ) = R(t w ) g(C) . (10) 
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Figure 5. The two-time correlation length, £, in the RFIM. (a) £ as a function of 
time-delay, t— t w for several values of t w given in the key at T = 1 and H = 1. (b) £ as a 
function of the global correlation in a parametric plot at T = 1 and H = 1 . (c) Scaling 
£{t,t w ) = R(t w ) g(C) at two temperatures and two values of the random held using 
three waiting-times t w for each set of parameters. The pure case H = 0, T = 1 is also 
included with a very short t w to avoid equilibration. 



The limit g(C = 1) = is found by taking t = t w , that corresponds to C — 1 [extending 
the scaling form ( ITUl) to include the thermal regime]. In this case C^(r;t,t) = 1. If 
one uses C 4 (r;t, t) = C 4 (r/£, C(t, t) = 1), see Sect. I3.1.4[ then must vanish 

to obtain C 4 independent from r, and this imposes g(l) = 0. In the other extreme, 
when t t w and C = one expects g(0) = 1. The reason is the following. 
lim t>4u) C 4 (r; t, t w ) = C^r, t)C2(r, t w ), for the temporal decoupling of C 4 can be done in 
the t 3> t w limit. Recalling that C^ir, t) oc f2{r/R(t)) with \\m x ^o f 2 (x) = 1, the only 
spatial contribution to lim t>tiu C 4 (r; t, t w ) comes from the term C^r, t w ) oc f2{r/R{t w )). 
Using \\m t ^,t w £,(t, t w ) = R(t w )g(0) and further assuming that the functional forms of 
C±{x) and fii^) are, to a first approximation, the same we deduce g(0) = 1. 

Figure [5] (c), where we plot £,(t,t w )/ R(t w ) versus 1 — C(t,t w ) for different t w , 
illustrates the validity of the scaling hypothesis (TTOl) . We see that, as expected, 
g(C = 1) = and it seems plausible that \im c ^o g(C) = 1. The scaling function g 
is found to satisfy super-universality, i. e. it is independent of H and T. 



3.1.4- C4 and super-universality. Using the monotonicity properties of C as a function 
of t — t w and t w , and of £ as a function of t w and 1 — C we can safely exchange the 
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Figure 6. (a) Test of scaling, Ci(r/^,C), and the super-universality of C4 for the 
parameters T and given in the key. Times t and t w are chosen in such a way that 
C(t, t w ) = 0.57 in all cases, (b) The same data in linear-log scale showing that C4 — C 2 
is very close to an exponential at short r/£. 



dependence of C4 on the two times by a dependence on £ and C. In other words, 
C^r, £,C) where, again for simplicity, we did not write explicitly the dependence on T 
and H. Now, a reasonable scaling assumption is that one can measure r in units of £ 
such that 

C 4 (r, t, t w ) = Gi(r/£(f , *«,), C(t, * ro )) . (11) 

In Fig. [6] we put this scaling form to the test and we examine the possible super- 
universality of C4. We use different values of the parameters t, t w , T, H such that 
C = 0.57 in all cases. Both scaling and super-universality relations are well satisfied. 
Note that the scaling relation in Eq. (fill) can also be transformed into 

C 4 (r; t, t w ) = C 4 (r/R(t w ), R(t)/R(t w )) (12) 

by using Eq. (jHJ). This last scaling form was also found for the O(N) ferromagnetic model 
in the large N limit although the scaling function does not have a simple exponential 
relaxation [20] . 



3.2. 3d EA 

A detailed analysis of the relaxation properties of similar correlations in the 3d EA model 
appeared in [TJ]. The spatial one-time correlation, C*2(r, t), vanishes identically in this 
model due to the quenched random interactions. The two-time self-correlation satisfies 
scaling with R ~ t 1 /^( T ) cannot be simply associated to a typical radius of equilibrated 
domains. The question as to whether the scaling function / is super-universal is not 
well posed since the T-dependent power l/z(T) can be absorbed in /. The four-point 
correlation allows for the definition of a two-time growing length scale £ that behaves 
qualitatively as in Eq. fflOl) . In Fig. [7] we present £(t, t w ) for the 3d EA. Its behaviour is 
very similar to the one of the RFIM exposed previously, but we would like to stress the 
fact that this quantity reaches much lower values in the 3d EA case (around 2a) than 
in the RFIM (around 15a). Figured (c) demonstrates that the superscaling property 
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Figure 7. Study of the two-time correlation length in the 3d EA model, (a) £ 
as a function of time-delay t — t w for several t w 's given in the key at T/T g = 0.6. 

(b) Evolution of £ with the global correlation in a parametric plot at T/T g — 0.6. 

(c) Test of the scaling hypothesis t w ) = i?(i lu ) g(C) with iJ(t) cx i 003 at T/T g = 0.3 
and T/T g = 0.6. 



does not hold in the 3<i EA model. We used R(t) oc t 03 for both temperatures and 
the resulting g(C) curves are significantly different. It is important to remark that no 
T-dependent power-law in R would make the two curves collapse. Turning back to the 
scaling of the two-time correlation and fixing the power law, C oc f[{t/t w ) om ] one finds 
f(x) ~ x~ 4 ' 5 (at T/Tg ~ 0.6) a much faster decaying power than in the RFIM. Note that 
previous estimates of the dynamic exponent using the one-time replica overlap [121 [13] 
yield l/z(T = 0.3T g ) w 0.045 a slightly larger value; the reason for the discrepancy 
could be traced to the lack of accuracy in the determination of £ and then R. 

3.3. Colloidal glasses 

The structure factor of colloidal suspensions and Lennard- Jones mixtures are obviously 
very different from the one of a sample undergoing ferromagnetic ordering. Still, two- 
time self-correlations satisfy scaling with R(t) cx t x l z although a clear interpretation of 
R is not available. 

Castillo and Parsaeian studied £ in a Lennard- Jones mixture of particles undergoing 
a glassy arrest. One notices that, at short time delays (t — t w ~ 10 molecular dynamic 
units), £ is monotonic with respect to t — t w and t w in this system, while one needs to 
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reach much longer time delays (and indeed go beyond the simulation window) in the 3d 
EA and RFIM cases [cfr. Figs. [5] (a) and [7] (a) to the first panel in Fig. 2 in [17]]. A 
form like ( jTUl) describes £ in this case too with R(t) ~ t 1/z and l/z ~ 0.1. 

The two-time correlation length of colloidal suspensions was analysed in [18] using 
a mapping to a spin problem. The data for £ remains, though, quite noisy and although 
a similar trend in time emerges the precise functional form is hard to extract. 

3.4- Summary 

In short, the macroscopic correlations in all these systems admit the same dynamic 
scaling analysis although there is no clear interpretation of R as a domain size in the 
case of the 3d EA and colloidal suspensions. 

4. Fluctuations 

An approach apt to describe problems with and without quenched randomness focuses 
on thermally induced fluctuations [19]. The local dynamics can then be examined by 
studying two-time spin-spin functions which, instead of being spatially averaged over 
the whole bulk, are only averaged over a coarse-graining cell with volume V r = (2Z) 3 
centered at some site r [16] : 



One can then characterize the fluctuations by studying their probability distribution 
function (pdf) p(C r ;t,t w ,l, L,T, H) with mean value C(t,t w ). 

In general, the variation of p(C r ) with the size of the coarse-graining boxes is as 
follows. For I < R the pdf is peaked around qea and has a fat tail towards small values 
of C r including negative ones. Indeed, well in the coarsening regime, most of the small 
coarse-grained cells fall inside domains and one then expects to find mostly a thermal 
equilibrium distribution - apart from the tail. For larger values of I such as I ~ R, a 
second peak close to C appears and the one at ^ea progressively diminishes in height. 
For still larger values of Z, the peak at qea disappears and a single peak centered at C 
(the mean value of the distribution) takes all the pdf weight. 

At fixed temperature and field, the pdf p(C r ;t,t w ,l, L) in the RFIM depends on 
four parameters, two times t and t w and two lengths I and L. In the aging regime the 
dependence on t and t w can be replaced by a dependence on C(t,t w ) and t w ), the 
former being the global correlation and the latter the two-time dependent correlation 
length. Indeed, C(t, t w ) is a monotonic function on the two times [cfr. Fig. [3] (a)] and £ 
is a growing function of t (cfr. Fig. [[]), thus allowing for the inversion (t,t w ) — > (C, £). 
Note that we do not need to enter the aging, coarsening regime to propose this form. 
One can now make the natural scaling assumption that the pdfs depend on £, the coarse- 
graining length /, and the system linear size L through the ratios //£ and l/L. In the 




(13) 
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Figure 8. Pdf of local two-time functions C r in the RFIM at T = 1 and H = 1. The 
waiting-times are given in the key and time t is chosen such that C(i, t w ) — 0.6. (a) C r 
is coarse-grained on boxes of linear size I = 9. (b) C r is coarse-grained on boxes with 
variable length I so as to keep l/£{t, t w ) ~ 0.7 constant. The collapse is much improved 
with respect to panel (a). 




Figure 9. Test of the scaling and super-scaling hypothesis. The two pairs of t and t 
are the same as in Fig. Hand C = 0.6 as well, (a) Z/£ ~ 1.4. (b) l/£ ~ 2.9. 



end, the pdfs characterizing the heterogeneous aging of the system read 

p(C r ;C(t,t w ),l/£(t,t w ),l/L) . (14) 

We numerically test this proposal by assuming that the thermodynamic limit applies 
and the last scaling ratio vanishes identically. Figure El (a) shows the pdfs at two pairs 
of times t and t w such that the global correlation C(t,t w ) is the same, and 1 = 9. It is 
clear that the two distributions are different. In panel (b) we further choose / so that 
l/£ ~ 0.7 is also fixed. The two distributions now collapse as expected from the scaling 
hypothesis Eq. (fl4l) . Note that another peak at C = — 1 exists, though with a lower 
weight. Figure M (a) and (b) show the scaling for //£ ~ 1.4 and l/£ ~ 2.9, respectively. 
While the collapse is still good in the case of panel (a), it is not satisfactory in panel (b). 
Indeed, this plot suffers from the fact that the thermodynamic limit is far from being 
reached (l/L ~ 0.15 is not so small). 
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Figure 10. Pdf of local correlations in the 3d EA model at T/T g = 0.3 and T/T g = 0.6, 
for two waiting-times t w such that C = 0.6 and //£ = 2.9. The solid line (red) displays 
the super-universal pdf in the RFIM. 



In Fig. O (a) we used several values of T and H and we found that all pdfs collapse 
on the same master curve. We conclude that as long as coarse-graining lengths are not 
too close to the system size, the pdf of local correlation satisfy the scaling (fill) with a 
scaling function that is super-universal. 

Let us now compare the forms of the pdfs in the RFIM and 3d EA model. In the 
RFIM the peak at oea is visible until //£ ~ 2. Given that in this model £ is quickly 
rather large (£ reaches 15a in the simulation time-window) one has a relatively large 
interval of I for which the peak at oea can be easily seen. Instead, in the 3d EA the 
two-time correlation length grows very slowly and reaches only £ ~ 2a in similar times, 
meaning that the peak at oea is hardly visible as soon as one coarse-grains the two-time 
observables [Ej. 

Figure flOl demonstrates that the pdf of local correlations is not super- universal with 
respect to T in the 3^ EA model, and compares the functional form at two temperatures, 
T/Tg = 0.3 and T/T g = 0.6, with the one in the RFIM. The global correlation, C, 
and the ratio of coarse-graining to correlation lengths, //£, are the same in all curves. 
Although qualitatively similar, the pdf in the RFIM and 3d EA models are different, 
with the RFIM one being more centered around the global value. 

The study of Lennard- Jones mixtures in [23] used a constant coarse-graining length 
and the pdfs of local correlations at constant C showed a slow drift that should be cured 
by taking into account the variation of £. In colloidal suspensions the scaling form (1131) 
is well satisfied [18]. In the context of coarsening phenomena these pdfs are to be 
compared to the ones calculated for the O(N) model in its large N limit [20] . 

5. Conclusion 

We performed an extensive analysis of the dynamics of the RFIM in its coarsening 
regime. We showed that the equal-time correlation functions, global two-time correlation 
functions, and the four point correlation functions obey scaling and super-universality 
relations in the aging regime. The scaling relations, by means of the typical growing 
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length, R oc lnt/r, reveal a non-trivial time-invariance for these statistical objects. 
Super-universality encodes the irrelevance of quenched randomness and temperature on 
the scaling functions and it is demonstrated by the fact that they are the same as for 
the pure Ising case. 

In the 3d EA, similar scaling forms were found for global two-time correlations 
and four-point correlations [H]. The function R(t) could be associated to a domain 
radius though a clearcut confirmation of this is lacking. On the contrary, the results of 
recent large scale simulations have been interpreted as evidence for an SK-like dynamic 
scenario [15]. The one-time function playing the role of the domain radius is a very 
weak power law, t a03 at T/T g ~ 0.3 — 0.6, and, in consequence, the two-time correlation 
length reaches much shorter values than in the RFIM in equivalent simulation times. 
Super-universality (with respect to temperature) does not apply in this case. 

A similar scenario applies to the Lennard- Jones mixtures [23] and colloidal 
suspensions [18]. The two-time correlation length remains also very short in accessible 
numerical and experimental times. 

In all these systems the analysis of local fluctuations of two-time functions leads to 
scaling of their probability distribution functions. In the RFIM these also verify super- 
scaling with respect to T and H. In the 3d EA they do not. The intriguing possibility 
of a kind of super-scaling in colloidal suspensions (with respect to concentration) has 
been signaled in [18J and deserves a more careful study. 

We conclude that all these systems, with a priori very different microscopic 
dynamic processes admit a similar dynamic scaling descriprion of their macroscopic 
and mesoscopic out of equilibrium evolution. 
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